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ABSTRACT 

We study the chemical abundances of the interstellar medium surrounding high z gamma-ray bursts 
(GRBs) through analysis of the damped Lya systems (DLAs) identified in afterglow spectra. These 
GRB-DLAs are characterized by large H I column densities A'hi and metallicities [M/H[ spanning 1/100 
to nearly solar, with median [M/H[ > —1 dex. The majority of GRB-DLAs have [M/H[ values exceeding 
the cosmic mean metallicity of atomic gas at z > 2, i.e. if anything, the GRB-DLAs are biased to larger 
metallicity. We also observe (i) large [Zn/Fe[ values (> -1-0.6 dex) and sub-solar Ti/Fe ratios which imply 
substantial differential depletion, (ii) large a /Ye ratios suggesting nucleosynthetic enrichment by massive 
stars, and (iii) low 0*^/0+ ratios (< 10"''). Quantitatively, the observed depletion levels and 0*^/0+ 
ratios of the gas are not characteristic of cold, dense H I clouds in the Galactic ISM. We argue that the 
GRB-DLA represents the ISM near the GRB but not gas directly local to the GRB (e.g. its molecular 
cloud or circumstellar material). We compare these observations with DLAs intervening background 
quasars (QSO-DLAs). The GRB-DLAs exhibit larger A'hi values, higher a/Fe and Zn/Fe ratios, and 
have higher metallicity than the QSO-DLAs. Although these differences are statistically significant, 
the offsets are relatively modest (Ahi excepted). We argue that the differences primarily result from 
galactocentric radius-dependent differences in the ISM: GRB-DLAs preferentially probe denser, more 
depleted, higher metallicity gas located in the inner few kpc whereas QSO-DLAs are more likely to 
intersect the less dense, less enriched, outer regions of the galaxy. Finally, we investigate whether dust 
obscuration may exclude GRB-DLA sightlines from QSO-DLA samples; we find that the majority of 
GRB-DLAs would be recovered which implies little observational bias against large TVhi systems. 
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iFiore et Hifm^ iChen et al. l[2005t iProchter et aTEoO^ 
and transitio ns from the ISM surrounding the GRB event 
(e.g. iSavaglio. Fall fc Fiord 120031: IVreeswiik efoL] 12004 
iProchaska et al. I 20061 ). And, analogous to absorption 
line studies of quasars, the se obse rvations reveal the H I 



1. INTRODUCTION 



For the past decade, high resolution spectroscopy of 
distant quasars (QSOs) have enabled detailed studies of 
the chemical abundances of high z galaxies (e.g. | Lu et al. 



1996: Proc haska fc W olfe 1999; D essauges-Z avadsk v et al. , ^^^^ 

2001; ,Molaro et al. 112001. : Prochaska et al. Ii2001; Tedoux. Petitiee^^^^Srm najgd 
2003[) . These properties of the interstellar medium (ISM) ' 



(Vrceswii k et al. 1200 4: Jakobsso n et al. I 



have been surveyed extensively via the damped Lya sys- 
tems (DLAs), galaxies ii itersected by sightlines to distant 
quasars (for a review, see lWolfe, Gawiser fc Prochaskal20()5l ). 
These QSO-DLA sightlines sample galaxies according to 
gas cross-section, i.e. the outer regions of a galaxy's ISM 
are more frequently probed because differential cross-section 
is proportional to radius. If star formation (SF) occurs in 
compact regions in high z galaxies (i.e. molecular clouds), 
then Q SO-DLAs would only rarely probe these environ- 
ments (IZwaan fc Procha ska 2006). This conclusion holds 
true independent of any additional biases rela ted to dust 
obscu ration by the SF regions themselves (e.g. I Fall fc Peil 
I1993D . QSO-DLA surveys are an inefficient approach to 
directly studying the gas in high z SF regions. 

Like quasars, gamma-ray bursts (GRBs) provide a bright 
- albeit transient - light beam from the distant universe. 
Imprinted on the roughly power-law spectrum of the after- 
glow are the signatures of the intergalactic medium (e.g. 



TTiretauicity fjFvnbo et al. 1 120061 : IProch aska 2006 
chemical abundances, differentia l depletion (ISavaglio fc Falll 



I2004D . and kinematics (Paper II: [Prochaska et al. II2007I ) of 
the gas along the sightline. 

The progenitors of long-duration GRB are believed to 
be massive stars residing in active sta r forming regions 
(lWooslevlll99"3t iWooslev fc bIoo^[2006I) . At low redshift, 
this association is well established: supernova h ave been 
identified at the sam e position a s GRB events ( e.g. lHiorth et al. I 
2003: Stanek et al. ^ 2003: Mirabal et al. 1120061 but see also 
Fvnbo et al. I (|2006D ). Furthermore, GRBs have b een found 



recen tly to be located within Wolf-Rayet galaxies (jHammer et al. I 
l2006f) . The connection between GRB and star-forming re- 
gions is inferred at high redshift. GRBs are found exclu- 
sively in galaxies that are blue and show nebular emission 
lines indicative of ongoing star formation (|Le Floc'h et al. I 
|2003|) . Furthermore, high precision astrometry of long- 
duration GRB afterglows reveal these events occur exclu- 
sively within 10 kpc of the flux-weighted centroid of their 



2 



host galaxy (jBloom. Kulkarni fc Diorgovskil2002HFruchter et g^d referentiallv probe the outer regions of galaxies whereas 



l2006t ). Therefore, the sighthnes to GRBs should probe the 
ISM within or at least near star-forming regions. In this 
respect, studies of the ISM along sightlines to GRBs com- 
plement those toward quasars. 

The first obse rvations of GRB-DLAs indicated large H I 
column density (jCastro et al. l2003tlVreeswiik et al. Il2004[ ) 
and la rge metal- line equivalent widths and co lumn den- 
sities (|Metzger et on 119971: ISavaglio. Fall fc Fiore 20033 . 
These data suggested modestly enriched gas (> 1/10 so- 
lar) and substantial depletion levels but the analysis is lim- 
ited by line-saturation (i.e. lo w-resolution spec troscopy of 
lines with large optical dept h: iPro chaska 2006) . With the 
launch of the Swift satellite (jGehrela, 2000 1 , GRBs are de- 
tected at a rate of « 2 per week with rapid, precise localiza- 
tions enabling cchcUe observations of a modest sample of 
GRB-DLAs ( Chen et al. 20 05; Proch aska, Ch en & Bloo^ 



2006|:IVreeswiik et I20Q6I: iProchaska et aZ. 120061: iPiranomonte et al 



GRB-DLAs are gene rally located within the inner few kpc 
of their host galaxies ([Bloom. Kulkarni fc Diorgovskil2002[ ) . 
This could imply higher A''hi values, metallicity, depletion 
levels, molecular fractions, and differential rotation along 
the GRB-DLAs sightlines. Therefore, a secondary theme 
of this paper is to test the hypothesis that GRB-DLAs and 
QSO-DLAs are drawn from the same parent population of 
high z galaxies and that any differences in ISM properties 
are explained by their average impact parameters through 
the galaxies. 

This paper is organized as follows. Section [2] presents 
the experimental design, and observational samples. In 
§ [31 we characterize the H I column densities, metallic- 
ities, relative abundances, and atomic carbon abundance 
of GRB-DLAs. We discuss the implications of these results 
land conclude with a list of future directions. 



20071) . These observations provide precise column density 



measurements of many metal-line transitions and permit 
the analysis of the chemical abundances in the gas sur- 
rounding GRB. 

In this paper, we describe the chemical abundances of a 
modest sample of damped Lya systems associated with 
the ISM surrounding gamma-ray bursts. We focus on 
echelle observations acquired in the past two years, but 
also include pre-Swift observations. Our principal goal 
is to describe the physical conditions of the interstellar 
medium within and/or near star formation regions in the 
young universe. To frame the discussion, we will compare 
the GRB-DLA observations against similar sets of obser- 
vations for damped Lya systems along quasar sightlines 
(QSO-DLAs). 

There are at least three reasons why one may expect 
the ISM characteristics of the GRB-DLAs to differ with 
those of QSO-DLAs. First, the two samples may arise 
from overlapping yet non-identical populations of high z 
galaxies. By definition, the QSO-DLAs correspond to 
galaxies with large H I surface densities. These are se- 
lected independent ly of any emission or stellar property 
([Wolfe et al. I [19861 ) and cosmological simualtions suggest 
a s ubset may not even be l ocated within dark matter ha- 
los ([Razoumov et al. [[2006() . Long-duration GRB, in con- 
trast, are known to occur only within star- forming galaxies 

and are expected to roughly trace the on-going star forma- , , j /-vor-. r^-.^ rnrxi r /-irm i j. 

J. /orn-r. I rn X - l li nnni rn ■ V, ■ " rr^ , i o r>i b sckgrouud QSO. 1 hc U V hux trom a GRB aiterglow that 

tionrate (SFR:e.g. Totamlll999HRamirez-Ruiz. Trentham fc Blami , • , r , • , .... ° 

r^i — ^ — 1^ Innnrj; — Tir 7 — 7^ — ■ ] IS" typical 01 spectroscopic observations (i.e. peak magni- 

2002t Guetta fc Piranl2007n . Of course, star-forming galax 



EXPERIMENTAL DESIGN AND OBSERVATIONAL 
SAMPLES 

2.1. Sightline Geometries 

Before discussing the observational data, it is impor- 
tant to comment further on a few fundamental differences 
between studying DLAs along GRB sightlines versus back- 
ground QSOs. The most obvious differences are that the 
GRB sightline originates within the galaxy giving rise to 
the damped Lya system and, almost certainly, from within 
a young star-forming region. This implies several impor- 
tant consequences for the DLAs probed by GRB after- 
glows. 

First, because the phenomena originates from within the 
DLAs themselves, a GRB sightline intersects less gas col- 
umn than a background QSO sightline would at the same 
position on the sky. Therefore, one underestimates the 
average H I column density and may underestimate the 
total velocity field that occurs along the sightline. Sec- 
ond, it is possible that gas local to the progenitor, i.e. 
circumstellar material and/or molecular cloud gas, con- 
tributes significantly to the GRB-DLAs (but see below). 
One may expect especially large dust-to-gas ratios, high 
molecular fractions, and/or peculiar chemical abundance 
ratios. Third, the GRB afterglow radiation can signifi- 
cantly affect the DLAs observed along its sightline whereas 
QSO-DLAs are selected to have large distance from their 



ies presumably also have large H I column densitieq^, as 
evidenced by the DLA profiles in GRB afterglows. 

Second, even if QSO-DLAs and GRB-DLAs are drawn 
from the exact same parent population of host galaxies, 
the QSO-DLAs are selected according to H I covering frac- 
tion on the sky while GRB are expected to track cur- 
rent SFR. One may expect QSO-DLA to preferentially 
arise in galaxies with extended H I disks (e.g. low-surface 
brightness) whereas GRB p referentialy occu r within high- 
surface brightness galaxies ([Mo. Mao fc W hite 1998). Fi- 
nally, and perhaps the most important, QSO-DLAs will 

^Similarly, we may expect that all galaxies with large H I col- 
umn density will exhibit star formation. Indeed, the presence of 
heavy elements in all Q SO-DLAs indicates previous star formation 
(|Prochaska et 0171 12003) and C II* fine-structure absorption implies 
ongoing star formation (Wolfe. Prochaska fc Gawiscr.2003,) . 



typical of spectroscopic observations (i.e. peak 
tude of Rk. 15 at early times) is sufficient to ionize an H I 
column of approximately 10^^ cm~^ a t 5 pc fro m the after- 
glow in only a few minutes time (e.g. Draine fc Hao 2002; 
[Perna fc Lazzatil [2002[ : [Prochaska. Chen fc Bloonj [20061 ). 
Similarly, theoretical treatments predict the X-ray and UV 
components of the afterglow can destroy dust and molecu- 
lar hydrogen out to 10 to lOOpc ( f Wax man fc Draind'2000^: 
[Fruchter. Krolik fc Rhoadi[200l[ i [Draine fc Hao 2002). If 
neutral gas and dust existed at distances of ~ lOpc prior 
to the burst, it is likely removed by the afterglow before 
spectroscopic observations are initiated. In this respect, 
the contribution of gas local to the GRB is somewhat min- 
imized by the GRB event itself. By the same token, gas 
local to the star forming regions may have been attenu- 
ated by the afterglow radiation field and our observations 
may not precisely reflect the conditions of this gas in the 
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Table 1 
GRB-DLA SAMPLE 



GRB 


RA 


DEC 


ZGRB 


Instrument 


R 


Mg I" 


Exc. Fell" 


Ref 


GRB990123 


15:25:30.34 


+44:45:59.1 


1.600 


Keck/LRIS 


1,000 


Y 


N 


1 


GRB000926 


17:04:09.00 


+51:47:10.0 


2.038 


Keck/ESI 


5,000 


Y 


N 


2 


GRB010222 


14:52:12.55 


+43:01:06.2 


1.477 


Keck/ESI 


5,000 


Y 


N 


3 


GRB011211 


11:15:17.98 


-21:56:56.2 


2.142 


VLT/FORS2 


1,000 


? 


N 


1,4 


GRB020813 


19:46:41.87 


-19:36:04.8 


1.255 


Keck/LRIS 


1,000 


Y 


Y 


5 


GRB030226 


11:33:04.93 


-f 25:53:55. 3 


1.987 


Keck/ESI 


5,000 


7 


N 


6 


GRB030323 


11:06:09.40 


-21:46:13.2 


3.372 


VLT/FORS2 


1,000 


Y 


N 


7 


GRB050401 


16:31:28.82 


+02:11:14.8 


2.899 


VLT/FORS2 


1,000 


7 


N 


8 


GRB050505 


09:27:03.20 


+30:16:21.5 


4.275 


Keck/LRIS 


1,000 


7 


N 


9 


GRB050730 


14:08:17.14 


-03:46:17.8 


3.969 


Magellan/MIKE 


30,000 


7 


Y 


10 


GRB050820 


22:29:38.11 


-1-19:33:37.1 


2.615 


Keck/HIRES 


30,000 


N 


N 


11 


GRB050904 


00:54:50.79 


+14:05:09.4 


6.296 


Subaru/FOCAS 


1,000 


7 


7 


12 


GRB050922C 


19:55:54.48 


-08:45:27.5 


2.199 


VLT/UVES 


30,000 


W 


Y 


13 


GRB051111 


00:08:17.14 


-00:46:17.8 


1.549 


Keck/HIRES 


30,000 


Y 


Y 


14,11 


GRB060206 


13:31:43.42 


-1-35:03:03.6 


4.048 


WHT/ISIS 


4,000 


7 


7 


15 


GRB060418 


15:45:42.40 


-03:38:22.80 


1.490 


Magellan/MIKE 


30,000 


Y 


Y 


11 



"Strong Mgl absorption (W=w eak). GRB-DLA with strong Mg I absorption are likely to have the majority of their gas at distances greater 
than 50pc from the afterglow llProchaska. Chen &: Bloomll2006t) 

''Positive detection of absorption from excited levels of Fe+. This measurement is difficult in low-resolution data and a non-detection should be 
interpreted as ambigous. The absence of absorption in higher resolution data, however, indicates the Fe+ gas (observed via resonance lines) is 
at large distance from the afterglow. 



References. — 1 
J2TO3); 6: Shin et a/. (2001 
[Prochaska et al. I II2006I') : 12' 
II2006I ) 



'Savaglio . Fall fc Fior d fSOOf); 2: 'Castro et al. ' ^200^; 3: 'Mirabal et~al. 
7: J^eg^wijk et al. 



IKawai et al. i I i2006ll : 



(2004); 8: Watson et al. (2006) 
13: iPiranomonte et al. i (i2007l ) 



||200^): 4: IVreeswiik et al I ll2006h: 5 : 'Barth e t aFI 
IBerger et al. I (120061); 10 : ICheneTaL (2005) ; 11: 
14: IProchaska. Chen fc Blooinl (120061) ; 15: \Fynho etaT] 



9: 



absence of a GRB. 

On the other hand, a side from the UV pumping of fine- 
struc t ure hnes dP rochas ka. Chen fc Bloomll2006l : iDessauges 
120061: IVreeswiik al. I I2006D. there is no compeUing case 
of a GRB afterglow having attenuated its surrounding 
medium. There have been no substantiated reports of hne 
variabihty in any resonance hne including, and most im- 
portantly, Mg I transitions. The majority of GRB-DLAs 
exhibit strong Mg I absorption coincident in velocity with 
the resonance and excited trans i tions of low-ion specie s 
(jProchaska. Chen fc BloM^l2006HProchaska et al. II2006D . 
This atom has an ionization potential IP=7.7 eV and would 
be ionized by the afterglow if it occur s with in « 100 pc of 
the event. IProchaska. Chen fc Blooinl (|2006f ) have argued, 
therefore, that the majority of neutral gas observed in 
GRB-DLAs is at a distance exceeding 100 pc from the pro- 
genitor. Meanwhile, the detection of UV pumped excited 
states of Si+, 0°, and Fe+ require that t he gas is located 
within ss 1 kpc of the afterglow. Indeed, IVreeswiik et al. I 
()2006^ have analyzed variations in the level populations of 
the excited levels of Fe"*" and Ni+ and constrain the dis- 
tance of the GRB-DLAs along the GRB 060418 sighthne 
to be at 1 ± 0.2 kpc from the afterglow. 

Table |T] indicates those GRB-DLAs that exhibit strong 
Mg I absorption. We also identify the GRB-DLAs which 
show absorption from excited states of Fe+. Those cases 
without Fe+ fin e-structure absorption must also l ie at dis- 
tance s > 100 pc (jProchaska. Chen fc Bloom.2006l : [Chen et 
[2003). These results indicate that the majority of gas in 
GRB-DLAs is not local to the star formation region and 
should not be adversely affected by the GRB afterglow. 
Our current expectation is that O and B stars from the 



Zavi 




Fig. 1. — This cartoon illustrates the likely differences between 
QSO-DLA and GRB-DLA sightlines. The former have randomly 
intersected a foreground galaxy. These QSO-DLA sightlines corre- 
spond to a cross-section selected sample and should preferentially 
intersect the outer regions of the ISM in high z galaxies. In con- 
trast, the GRB-DLAs are constrained to originate from within the 
ISM of their host galaxies, presumably the H II region produced by 
massive stars in a star-forming region. These GRB-DLA sightlines 
are expected (and observed) to originate within the inner few kpc 
of the ISM. 



star-forming region have ionized all gas within w 100 pc of 
the GRB producing large H II regions. Figure [T] presents a 
cartoon illustration of a GRB and a QSO sightline through 
a galaxy comprised of a neutral ISM surrounded by pre- 
dominantly ionized halo gas clouds. The GRB sightline 
originates near the center of the galaxy, within an H II re- 
gion associated with ongoing star- formation. In contrast. 
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Table 2 
QSO-DLA SAMPLES 



Sample 


Description 


Ndla 


Instrument (s) 


Ref 


HR-A 


High-resolution full sample (all) 


153 


HIRES, ESI, UVES 


1,2,3,4,5 


HR-S 


High-resolution 'statistical' sample 


112 


HIRES, ESI, UVES 


5 


HR-E 


High precision echelle 


71 


HIRES, UVES 


1,2 


SDSS 


True (A''hi) statistical sample 


475 


SDSS 


6,7 



References. — 1: IProchaska et a l. (2007); 2: Dessauges- Zavadskv et al. I (|200^): 3: ILedoux et al \ HwM ): 4: IHerbert-Fort et al. I |[20M) ; 5: 
IProchaska et al. I ll200a) : 6: IProchaska fc Herbert- Forti (.20041) : 7: IProchaska. Herbert-Fort fc Wolfd II2005I) 



the QSO-DLA sightlines will preferentially penetrate the 
outer regions of the ISM. 

Independent of the differences one may expect for QSO- 
DLAs and GRB-DLAs, there may be common trends (e.g. 
metallicity vs. dust depletion) characteristic of the ISM 
of all high z galaxies. By considering the QSO-DLAs and 
GRB-DLAs together, we hope to reveal these global trends 
in addition to highlighting differences in the samples. 

2.2. GRB-DLA and QSO-DLA Samples 

The QSO-DLA samples are drawn from two sources: 
(i) high-resolution echelle and ech ellette observatio ns of 
quas ars acquired with th e HIRES (jVogt et al. Ill994l) and 
ESI ([Sheinis et al. Il200dl ) spectrometers at the Keck Ob- 
servatory and the UVES (jPekker et al. I|2000D spectrom- 
eter at the VLT Observatory; and (ii) the low-resolution 
QSO-DLA surveys of the Sloan Digital Sky Survey (SDSS 



assuming solar metallicitj^; (ii) spectra with sufficient signal- 
to-noise and resolution to study the chemical abundances 
of the gas. The sample is also limited to GRB-DLAs where 
we could access the data or where precise equivalent width 
measurements were reported in the literature. Table[l]lists 
the GRB-DLAs comprising our sample and describes the 

spectral observations. 

As with QSO-DLAs (|Fall fc Peil 119931 ) . it is hkely that 
dust obscuration plays a role in defining the GRB-DLA 
sample. There are several examples in the literature of 
highly reddened afterglows which are likely extinguished 
by dust in the GRB host galaxy (e.g. iLevan et al. 1 120061 : 
iPellizza et al. ifJoOQ ). At present, no comprehensive study 



on the incidence of optically 'dark' bursts has been per- 
formed nor an evaluation of the fraction of dark bursts 
which are cases of dust obscuration (as opposed to high 



z events). We expect, however, that dust obscuration is 
Prochaska fc Herbert-Fortl200l: IProchaska. Herbert-Fort fc Wo^W P"^^^'^^ ^° defining GRB-DLA samples. 

2.3. Absorption-Line Metallicities 

Metallicity is the mass density in metals relative to hy- 
drogen and helium gas. It is likely that oxygen dominates 
this quantity in nearly every astrophysical environment, 
with carbon and nitrogen secondary. Although these three 
elements exhibit low-ion transitions, the majority are ei- 
ther too strong, too weak, and/or located within the Lya 
forest. Therefore, the metal abundance of the interstel- 



2005h . We will restrict the samples to zdla > 1.6 and 
DLAs which are greater than 3000 km s~^ from their back- 
ground quasar. The high-re solution sample (HR-A) i s 
summarized in these pa pers: Herbert-Fort et al. (2006); 
Ledoux et al. I (1 2006): D essauges-Zavadskv et al. I ()2006[ ): 
Prochaska et al. i (|2007|) . IProchaska et al.' { 2007f) have em- 



phasized that the full dataset is a heterogeneous sample of 
QSO-DLAs including a nu mber of systems selected o n the 
basis of strong metal-lines (jHerbert-Fort et al. II2006I) or as 

prom ising candidates for H':. absorption (|Ledoux. Petitiean fc Srlananil iium is frequently gauged by other elements. Un 



120031 ). As such, the H I distribution /(TVhi) of the high 
resolution sample does not follow the statistical distribu- 
tion derived from the SDSS. 

For a number of the comparisons made in this paper, 
we will take subsets of these samples. In particular, we 
define a pseudo-statistical sample of metallicity measure- 
ments (HR-S) w h ich is restricted to the compilation of 
IProchaska et al. I (|2003l ). We also construct a sample of 
high quality echelle-only observations (HR-E) for analysis 
related to relative abundances (e.g. a/Fe). The echellette 
observations, in particular, have too poor data quality to 



fortunately, many of the transitions which are accessible 
to analysis arise from elements that are refractory. To 
avoid the complications of depletion corrections (which 
can be an order of magnitude or more), observers have 
focused on non-refractory or mildly refractory elements. 
These include S, Si, and Zn. The latter is a trace ele- 
ment, log(Zn /H)(:) = —7.3, with a n uncertain nucleosyn- 
thetic origin ( Hoffman et al. \\l99^ and, therefore, should 
be considered cautiously. 

A point of great interest in GRB studies is the metallic- 
ity of the progenitor especially in the context of the col- 



gener allv achieve better than 0.1 dex precision (Prochaska et al. I lapsar model (Wooslcy 1993). In the coUapsar paradigm, 



I2OOI and are not included. Although the DLAs in these 
high-resolution subsets do not follow /(iVni) for a random 
sample, they were selected only on the basis of a large 
H I column density. Table [2] summarizes the various QSO- 
DLA samples considered in this paper. 

Our selection criteria for the GRB-DLA sample are (i) 
the presence of a damped Lya system (A^hi > 2 x 10^° cm~^) 



or a low- ion column density that requires TVhi exceed 2 x 10 ° cn ED65. : .Prochaska et al. 1,2005.) 



large angular momentum is required to power the GRB. 
Because high metallicity stars are expected to have signif- 
icant ma ss- loss by winds promot i ng the loss of angular mo- 
mentum (IVink fc de Ko teiJ2005!).[L anger fc NormanI (|2006t) 
and lWooslev fc Heged 112006) have argued that GRB pro- 

•3 We note that the few GRB sightlines with A^hi < 10^° cm~^ 
have very low column densities for low-ion transitions (IFiore et al. I 
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QSO-DLA 
GRB-DLA 



genitors will have low metallicity, i.e. less than 1/10 solar 
abundance. One may consider this to be a natural predic- 
tion of the coUapsar model. 

There is at least circumstantial evidence in support of a 
metallicity 'bias' from studies of z < 1 GRB host galaxies. 
The metallicities of the H II regions in a s mall sample of 
z < .5 galaxies hosting GR Bs are sub-solar (jProchaska et al. I 
[20P4';'SoliCT man et al. 1200 51 and several authors have noted 
that the values are systematically lower than the predicted 

distributi on for galaxies drawn randomly according to cur- 

rent S FR (jStanek et al. ll2006HKewlev et al. ll2007HModiaz et al. I 

l2007t ). Furthermore, the GRB host galaxies at z 1 are 

gener a lly sub-L^. consis t ent w ith low metallicity ijLe Floc'h et al. I 
l2003f ). iFruchter et"^ (|2006D compared a sample of z < 1 
GRB host galaxies against the host population for core- 
collapse supernova and demonstrated that GRB host galax- 
ies have systematically lower luminosity. The difference is 
approximately one magnitude which they suggested could 

be due to a metallicity bias. However, galaxies follow a „ t,, / , , ,x , , ■ ,■ ,r 

, IT , ■ ., 1 ;• ;i ; • 1 1 T 1 — 1 — W'^'ir ^r^=~Blue (dark, solid) curve traces the histogram oi A/jji 

metalhcity-lummosity relation that is roughly Imear (|Kobulnicky^^^K(>w.lc¥e GRB-DLAs. In comparison, we show a histogram 

of Nm values for QSO-DLAs drawn randomly toward background 
quasars (Prochaska. Herbert-Fort & Wolfe 2005) and normalized for 
presentation to have the same number of systems as the GRB-DLA 
sample. The GRB-DLAs have median value logA^ni = 21.7, which 
exceeds all but a few QSO-DLAs observed to date. The dotted 
red line traces the predicted distribution of A'hi values assuming a 
sample of sightlines originating inside an H II region located at the 
center of an H I exponential disk. The model shown has central 
column density logA^o = 22, scale-height to scale-radius h/R = 0.1, 
and an H II region with radius run = 4/i. 
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l2004f ) and a one magnitude difference in luminosity im- 
plies only an approximately 0.3 dex o ffset in metallicity. 
Recently, IWolf fc Podsiadlowskil ()2006f ) performed a thor- 
ough analysis showing the observations suggest a metallic- 
ity 'ceiling' for GRB progenitors but that this cutoff occurs 
at no lower than 1/2 solar. Therefore, while the offset in 
luminosities between the host galaxies of GRB and core- 
collapse supernova may be explained by a metallicity bias, 
this does not imply as severe an effect as promoted by 
studies of mass-loss. 

At high redshift, it is difficult to observe nebular lines to 
measure metallicities and our knowledge of galaxies is too 
poor to make robust statements from luminosity distribu- 
tions alone. Instead, one can measure the metallicity of 
the interstellar medium from absorption-line spectroscopy 
by comparing the total hydrogen column density with the 
total column density of metals. 

3. GAS PHASE ABUNDANCES IN GRB-DLAS 

In this section, we present column densities of hydro- 
gen and metals observed in the gas-phase in GRB-DLAs. 
These observations give the surface density and metallicity 
of the ISM, and the relative abundances reflect differential 
depletion and the underlying nucleosynthetic patterns in 
the gas. Finally, we comment on the abundance of atomic 
carbon. In all of this analysis, we draw comparisons be- 
tween the GRB-DLAs and QSO-DLAs. 

3.1. Hydrogen Column Densities 

In principle, the absorption-line spectra of rest-frame 
UV transitions (A < 3000 A) will give measurements of the 
atomic hydrogen and molecular hydrogen column densi- 
ties, iVni and N^i^. For a damped Lya system, the for- 
mer is best determined through a Voigt profile fit to the 
Lya A1215 transition. At A^hi > 10^° cm^^, the damping 
wings of the Lorentzian line-profile are well-resolved with a 
moderate resolution spectrum (FWHM< 5A). Therefore, 
aside from the redshift, the A^hi value is the most readily 
measured physical characteristic of a damped Lya system. 

Figure [5] presents the sample of A^hi measurements for 
o ur GRB-DLA sample su pplemented by the compilation 
of iJakobsson et oTI (12006 1) and compared against the dis- 
tribution of A^Hi values for QSO-DLAs measured from the 



Sloan Digital Sky Survey (iProchaska. Herbert- Fort fc Wolfd 

2005 ) . As previously reported (IVree swijk et al. 1120041 :1 IJakobsson et 
2006() . the GRB-DLA distribution is skewed to signifi- 
cantly higher A^hi values than QSO-DLAs. The difference, 
at least qualitatively, supports the hypothesis that QSO- 
DLAs probe the outer regions of galaxies whereas GRB- 
DLAs probe the inner, star-forming regions. 'Jako bsson et al~\ 
(j2006) discussed that the GRB-DLA distribution is roughly 
consistent with th at predicted for GRB em bedded within 
molecular clouds (jReichart fc Pric3l2002D . The authors 
note, however, that GRB exhibit too many sightlines with 
A'hi < 10^^ cm~^ compared to the prediction for molecu- 
lar clouds. They argue this may result from photoioniza- 
tion and /or because GRB are preferentially located at the 
edge (IWhite et aL \ 1999t) or even outside molecular clouds 
(jHammer et al. II2006D . 
As we discussed in ij |2.1i however, the presence of strong 



Mg I absorption in nearly every GRB sightline argues 
that the majority of GRB-DLA gas is located beyond » 
lOOpc of the event. This distance is comparable to only 
the lar gest giant molecul ar cloud complexes in the Local 
Group (|Blitz et al. II2007I) and it is reasonable to assume 
that the GRB-DLA gas is generally not associated with 
the molecular cloud hosting the GRB (see § 14. ip . Ex- 
amining Figure [5] in this light, the most salient question 
becomes: why do GRB-DLAs exhibit a preponderance of 
A'hi > 10^^ cm~^ measurements? Indeed, the gas mass 
required to average A^hi = 10^^ cm~^ at 100 pc along ran- 
dom sightlines is Mhi — IQ^ Mq; this exceeds the masses 
of even the la rgest molecular clouds in the Milky Way 
(jSolomon fc R ivolo 1989; Blitz 1993). It seems unhkely, 
therefore, that the observed H I gas corresponds to the 



Table 3 

GRB-DLA ABUNDANCE SUMMARY 



GRB 


Ic 




fa 
a 


[a/H] 




tc 


[Zn/H] 






[M/H] 




fe 


[Fe/H] 




[Ti/H]/ 


fc 


[N/H] 


a(N) 


iV(C°)^ 


GRB990123 




22" 













2 


-0, 


,97 


L.L. 


12 


-0 


.97 


L.L. 


2 


-2 


.00 


L.L. 













GRB000926 


21, 


.301°: 


,25 
,26 


2 


-1, 


.58 


L.L. 


1 


-0, 


,17 


0.15 


2 


-0 


.17 


0.29 


5 


-1, 


.49 


0.08 


-0.43 
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GRB010222 




22" 


2 


-1, 


.61 


L.L. 


2 


-1, 


,30 


L.L. 


12 


-1, 


.30 


L.L. 


1 


-2 


.01 


0.08 


-2.35 











GRB011211 


20, 


.401°: 


,20 
,20 


2 


-1, 


.36 


L.L. 











11 


-1, 


.36 


L.L. 


2 


-1, 


.80 


L.L. 













GRB020813 




22" 


2 


-1 


.31 


L.L. 


2 


-1, 


,17 


L.L. 


12 


-1, 


.17 


L.L. 


2 


-2 


.10 


L.L. 


-1.94 











GRB030226 


20, 


,501°: 


,30 
,30 


2 


-1, 


.31 


L.L. 


3 


-0, 


,47 


U.L. 


11 


--1 


.31 


L.L. 


1 


-1 


.05 


0.18 













GRB030323 


21, 


,901°: 


,07 
,07 


2 


-1, 


.56 


L.L. 


2 


-0, 


,87 


L.L. 


12 


-0, 


.87 


L.L. 


1 


-2, 


.40 


0.43 













GRB050401 


22, 


,601°: 


,30 
,30 


2 


-2, 


.16 


L.L. 


2 


-1, 


,57 


L.L. 


12 


-1, 


.57 


L.L. 


2 


-2 


.30 


L.L. 













GRB050505 


22, 


.051°: 


,10 
,10 


-4 


-1, 


.25 


L.L. 











11 


-1, 


.25 


L.L. 


2 


-1, 


.35 


L.L. 













GRB050730 


22, 


islg: 


,10 
,10 


4 


-2, 


.26 


0.10 











4 


-2 


.26 


0.14 


1 


-2 


.50 


0.12 




1 


-3.16 


0.10 


13.2 


GRB050820 


21, 


,001°: 


10 

,10 


4 


-0, 


.63 


0.04 


1 


-0, 


,71 


0.02 


4 


-0 


.63 


0.11 


1 


-1, 


.60 


0.09 


-0.80 


2 


-1.35 


L.L. 


12.7 


GRB050904 


21, 


,301°: 


,20 
,20 


-4 


-1, 


.10 


L.L. 











11 


-1 


.10 


L.L. 






















GRB050922C 


21, 


,601°: 


10 

,10 


4 


-2, 


.03 


0.10 


3 


-2, 


,02 


U.L. 


4 


-2 


.03 


0.14 


1 


-2 


.63 


0.01 


-1.60 


3 


-4.09 


U.L. 


12.9 


GRB051111 




22" 


2 


-1, 


.42 


L.L. 


2 


-0, 


,96 


L.L. 


12 


-0, 


.96 


L.L. 


1 


-2, 


.14 


0.01 


-2.23 









13.0 


GRB060206 


20, 


.851°: 


,10 
,10 


4 


-0, 


.85 


0.15 











4 


-0 


.85 


0.18 






















GRB060418 




22" 


2 


-1, 


.67 


L.L. 


1 


-1, 


,65 


0.04 


2 


-1, 


.65 


1.00 


1 


-2 


.24 


0.03 


-2.23 
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circumstellar medium which hosted the GRB. Instead, we 
contend the gas is associated with the nearby ISM of the 
galaxy. This hypothesis, however, must account for the 
large iVm values while allowing for an evacuated volume 
(i.e. H II region) with radius 100 pc or more. 

To explore this point further, we considered an H I disk 
described by a double exponential, 

"<^-«' = x=-(nr)=K^) • w 

with scale-height h, disk length Rd, and central H I column 
density Nq. We searched a wide parameter space of h/Rd, 
No, and radius of the H II region rHii- Over-plotted on 
Figure [5] is the prediction for random sightlines originating 
at the center and mid-plane (R — Z — 0) of an exponen- 
tial disk with logA^o = 22, h/Rd = 0.1 and rnii = 4/i. 
This is a reasonably good description of the observed dis- 
tribution {Pks = 11%); an ambient ISM characterized by 
logiVo = 22 with a large H II region {rnii > 2/i) produces 
a reasonable model for the observed A^hi distribution. It 
would be worthwhile to consider this model within the con- 
text of star forming galaxies in cosmological simulations. 

The association of GRB with star-forming regions raises 
the possibility that a significant fraction of the gas in GRB- 
DLAs is molecular. Conveniently, H2 gives rise to several 
band heads at A « 1000 A and one can directly measure 
the H2 column density. For optical spectroscopy, this re- 
quires ZGRB > 2 and blue wavelength coverage. It is also 
important to have relatively high resolut ion to disentangle 
the ab sorption lines from the Lya forest. iTumlinson et aU] 
()2007t ) have recently presented an analysis of five GRB 
sightlines and set an upper limit to the molecular fraction 
/(H2) < 10~^ in four of the systems. The only possi- 
ble detection is along the sightline to GRB 060206 with 
/(H2) w 10^3.5 (|Fvnbo et a/ . I I2006D. and even this valu e 
is likely to be an upper limit (|Tiimlinson" et g/. I[200l . 
Therefore, we will assume that A^H2 is a small fraction 
of the total hydrogen column density. These results lend 
further support to the interpretation of GRB-DLAs as be- 
ing dominated by ambient ISM instead of gas local to the 
progenitor. 

Lastly, there is the contribution to A^(II) from 11+ . Out 
to distances of several 10 pc, the star formation region (or 
GRB afterglow) will have produced an H II region with 
a potentially large A^(II+) value. That is, unless stel- 
lar winds from the GRB progenitor have nearly evacu- 
ated the region altogether. There is, of course, no di- 
rect means of measuring the H+ column density. Instead, 
we restrict our analysis of the metal column density to 
the atoms and ions that are dominant in neutral regions. 
These are termed "low- ions" - Si"*", Fe+, C+, 0°. We 
caution th at low-ions are not entir ely absent from H II re- 
gions fe.g. lHowk fc Sembachlll999l ). although the expected 
contribution would be small compared to neutral gas for 
A'hi > 10^^ cm^^. A detailed treatment of the ionization 
state of GRB-DLAs will be presented in a future paper. 
Here, we will assume that ionization corrections are small 
and comment on conclusions which are sensitive to this 
assumption. 

3.2. GRB-DLA Metallicities 

'^The following results are not too sensitive to the functional form 
of the radial profile. 



Ideally, the metal column density is derived from a sin- 
gle or set of unsaturated, resolved transitions of a low-ion. 
We will not attempt to measure the total column den- 
sity of each element (i.e. by summing all of the ionization 
states of a given element) but only the state which is dom- 
inant in neutral hydrogen regions (the low-ion). This can 
generally be achieved from spectra with high-resolution 
observations and modest signal-to-noise ratio. To date, 
however, many GRB spectra have been acquired with low- 
resolution spectrometers. These observations may provide 
accurate measurements of the equivalent widths (EWs) of 
strong lines, but precise column dens ities are diffi cult to 
derive due to modeling the line-profile (|Jenkinslll986f ). The 
high-resolution observations of GRB-DLAs, however, indi- 
cate that the metal-line profiles are comprised of multiple 
components (termed 'clouds') which exhibit a bimodal dis- 
tribution of column densities (Prochaska 2006). In most 
cases, the total column density is dominated by a single or 
few clouds whereas the equivalent width is the net sum of 
many weak clouds (Paper II). When the cloud column den- 
sities exhibit a bimodal distribution, the single-component 
COG analysis systematically underestimates the ionic col- 
umn densities (|Jenkinslll986l : iProcha ska 2006). In the fol- 
lowing, we will use low-resolution EW observations only 
to set lower limits to column density measurements and 
exclude these from relative chemical abundances. 

Table [3] summarizes the abundance measurements for 
the GRB-DLAs. We have corrected ions which exhibit 
fine-structure splitting of the ground-state (e.g. Fe+) by 
the observed contribution of these levels. In general, this 
implies an increment of < 0.05 dex. In five cases the H I 
column density and a non-refractory or mildly refractory 
metal abundance have been precisely measured. For the 
remainder of cases, we report a lower limit to the metal- 
licity because the metal abundance is a lower limit and/or 
there is no spectral coverage of the Lya transition. In the 
latter cases, we have assumed A^hi = 10^^ cm~^ which is 
approximately the mean logA'ni value of GRB-DLAs. 

We have evaluated the metallicity of the gas by adopting 
(in order of preference) the [S/H]3 value, the [Si/H] value, 
the [Zn/H] value, or the maximum of these if each has only 
a lower limit value. We prefer S and Si to Zn because these 
elements are w 1000 x more abundan t and because Zn ha s 
an uncertain nucleosynthetic origin (Hoffman et~dlil9M} . 
Again, we avoid highly refractory elements like Fe and Ni 
(even though they may dominate the opacity in stellar at- 
mospheres and therefore are particularly relevant within 
the collapsar model) because these elements may be de- 
pleted from the gas-phase. 

Figure [3] presents the metallicity measurements for the 
GRB-DLAs as a function of the age of the universe (set 
by zqhb) assuming the concordance cosmology (£7 a = 
0.72, rj„ = 0.28, /i = 0.73; Spcrgcl et a/. 2006). Over- 
plotted on the figure are the metallicity measurements for 
the statistical sample (HR-S) of QSO-DLAs. We restrict 
the discussion in this section to t < 3.5 Gyr (i.e. z > 1.65) 
where all of the GRB-DLAs have measured iVni values. It 
is apparent that both the GRB-DLAs and QSO-DLAs ex- 
hibit a large dispersion of values. The peak-to-peak range 
is « 2 dex, much larger than observational uncertainty. 
Second, both samples exhibit a metallicity 'floor' at ap- 

»[X/Y] = log(X/y) - \og{X/Y)Q 
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Fig. 3. — Metallicity [M/H] measurements for the GRB-DLAs (dark blue) as a function of the ag e of the universe corresponding to the 
observed absorption redshift and assuming the current concordance cosmology l lSpergel et al. ||2006I) . At z < 1.6 where Lya is lost below 
the atmosphere, the small arrows indicate lower limits to the metallicity assuming A^hi = lO^^cm"'^. At 2 > 1.6, the lower limits to 
[M/H] for the GRB-DL As are due to line-saturation. The lighter points show measurements for the QSO-DLAs assuming the HR-S sample 
i|Prochaska et al. Il2003i ^. The plus signs indicate A'^hi < 10^^ cm and the diamonds correspond to A^hi > lO'^^ cm ^. We also present the 
cosmic mean metallicity derived from the QSO-DLAs by taking the H I-weighted mean of the individual data points. Comparing the two 
distributions, we note that the majority of GRB-DLA values lie above the cosmic mean and that a significant fraction have [M/H] > —1. 



proximately 1/1000 solar abundance. This lower limit to 
the metallicities is also not observational; the sensitivity 
limit of the data is at least an order of magnitude lower. It 
remains an open question whether this floor is associated 
with early (PopIII) enrichment or rapid local enrichment 
in a ll galaxies exhibiting DLAs fWasserburg fc Qianl l200dl : 
iPro chaska et al. 2003; Qian & Wasserburg 2003). 

Also over-plotted in Figure [3] is the cosmic mean metal- 
licity of atomic gas. This quantity is derived from the 
AnT-weighted mean of QSO- DLA metallicities (jLanzettal 
119931 : IProchaska firZ1l20n3f ). Six of the 10 GRB-DLAs 
with Ahi measurements exceed the cosmic mean and sev- 
eral lower limits lie at the cosmic mean. Therefore, many 
and likely most of the GRB-DLAs have metallicities ex- 
ceeding the cosmic mean in the ambient ISM of high z 
galaxies. On these grounds, at least, the GRB-DLAs at 
high z do not appear to show a significant metallicity bias 
toward low values. 

Figure [3] presents a histogram of [M/H] values for the 
GRB-DLAs and QSO-DLAs (restricted to z > 1.6). The 
solid bars indicate GRB-DLA values while the dark open 



bars show lower hmits to [M/H] for the GRB-DLAs. We 
observe that that the GRB-DLA values roughly overlap 
the QSO-DLA distribution. A proper treatment of the 
observations is to perform a two-sample survival analysis 
(jFeigelson k NelsonI ITOSSI ) . Using the ASURV software 
package, we ran standard Gehan, logrank, and Peto-Peto 
tests using the Kaplan-Meier estimator and found the null 
hypothesis is ruled out at 99%c.l. We conclude that the 
metallicities of GRB-DLAs are larger than those for QSO- 
DLAs at z > 1.6. 

We note that many of the GRB-DLAs have a metallic- 
ity greater than 1/10 solar. At first glance, this appears to 
contradict assertio ns that GRB progenitors have metallic- 
itxless than 0.1 Zq (jLanger fc Normanll2006l : rWooslev fc Hegeil 
200 (1) . We will demonstrate in the next section, however, 
that the gas-phase abundances of S, Si, and Zn are all 
enhanced relative to Fe. The data allow, but do not re- 
quire, that the abundance of Fe is « 0.5 dex lower than 
the [M/H] values. It is possible, therefore, that at high 
z the GRB-DLA metallicities are not biased low and also 
that [Fe/H] < — 1 in most cases. Finally, the gas that we 
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Fig. 4. — Histogram of [M/H] measurements for the QSO-DLAs 
(light green) and the GRB-DLAs (dark) with both samples re- 
stricted to 2: > 1.6. For the GRB-DLAs, the open histogram traces 
the lower limits to [M/H] because of line-saturation. If these values 
are evaluated as measurements, then the two distributions are con- 
sistent with being drawn from the same parent population. If we 
increment the lower limits by -1-0.3 dex, however, the null hypoth- 
esis is ruled out at > 99%c.l. by a two-sided KS test. Similarly, 
a two-sample survival analysis l IFeigelson &: Nelsonlll985l) rules out 
the null hypothesis at greater than 99%c.l. 



observe is many lOpc's away from the progenitor and need 
not accurately reflect its metallicity. At present, however, 
we have no reason to suspect that the ISIVI values would 
systematically overestimate the GRB progenitor metallic- 
ity 

3.3. Relative Abundances 

We now turn our attention to the relative abundances 
of the GRB-DLAs. These are measured by comparing the 
gas-phase column densities of pairs of low- ions (X^, Y^) un- 
der the assumption that ionization corrections are small, 
i.e. [X/Y] = logTVx, -logA^r, -log(X/y)Q. With high 
resolution observations and moderate S/N data, one can 
frequently achieve 0.05 dex precision or better for logNx^- 
Because the observed ratios represent gas-phase abundances, 
however, the values reflect a combination of the underly- 
ing nucleosynthetic pattern and the effects of differential 
depletion onto dust grains. It is unfortunate that there 
is no element in the Fe peak which is non-refractory. As 
such, observers frequently employ Zn (a neighbor of the 
Fe peak) as a surrogate for Fe because it is non-refractory 
and because it traces Fe in the G alaxy at [Fe/H] > —2 
(jSneden. Gratton fc Crockeil ll99lD . We note, however, 
that Zn has an uncertain nucleosynthetic origin and should 
not directly trace Fe in galaxies whose star formation his- 
tory d iffers from the JVIilky Way (|Fenner. Prochaska fc GibsonI 
12004 ). Therefore, we consider it to be a surrogate for Fe 
but with a large systematic uncertainty. 

Throughout this section we restrict the QSO-DLA sam- 
ple to the high precision echelle measurements (sample 
HR-E; Table [ID. 



A key nucleosynthetic diagnostic of stars and galaxies is 
the a/Fe ratio, where a refers to the sequence of He fusion 
in massive stars, i.e. O, IVIg, Si, S, Ar. Standard theory of 
nucleosynthesis predicts these e lements are predominantl y 
produced by massive stars fe.g. IWooslev fc Weaver|[T995t ). 
Therefore, a comparison of the a abundance with the Fe- 
peak observations, whose production is dominated by the 
Type la SN of less massive stars, constrain the age and 
star formation history of the galaxy when plotted against 
gas metallicity (|Tinslevlll979D . Unfortunately, O and Mg 
are difflcult to measure via UV absorption lines because 
the transitions are either too strong or too weak. Instead, 
one typically estimates the a abundance with Si or S. In 
the QSO-DLAs, [S/Si] » and the eleme nts are roughly 
interchangeable (|Prochaska fc Wolfell200^ . For the GRB- 
DLAs, we have also adopted these two elements as the 
reference for a. The Fe-peak abundance, meanwhile, is 
determined from Fe, Ni, or Cr. 

Figure[5^ presents the observed a/Fe ratios for the GRB- 
DLAs against the QSO-DLAs as a function of a abun- 
dance, [a/H]. We reemphasize that these gas-phase abun- 
dances have not been corrected for differential depletion. 
Consider, first, the values for the QSO-DLAs. At low 
metallicity ([a/H] < -1.5), the QSO- DLAs follow a well- 
defined 'plateau' at [a/Fe]w 0.25 dex (jProchaska fc Wolfd 
[2002). This platea u matches th e one observed for metal- 
poor Galactic stars (jlVIc Williamii 1997) sug gesting the a/Fe 
enhancement has a nucleosynthet ic origin ( Lu et a/."1ll996l : 
iDessauges-Zavadskv et aZ. 1120061 ). At higher metalhcity, 



3.3.1. a/Fe 



the mean and dispersion of [a/Fe] rises, in stark contrast 
to the a/Fe trend in the IVIilky Way and that expected for 
other star formation histories (e.g. ISmecker-Hane et al. I 
l2002f ). The increase in a/Fe with increa sing metallicity is a 
clear signature for differential depletion ('Prochask a fc Wolfd 
[2002 ). i.e., large [a/Fe] ratios at high metallicity result 
from the greater adsorption of Fe onto dust grains. 

Turning to the GRB-DLAs, all of the [a/Fe] measure- 
ments are consistent with at least -1-0.5 dex and the ma- 
jority lie at greater than -1-0.6 dex. Two-sample survival 
analysis tests (e.g. Gehan) report a less than 0.1% proba- 
bility that the [a/Fe] values of the QSO-DLAs and GRB- 
DLAs are drawn from the same parent population. We 
conclude that the GRB-DLAs have systematically higher 
[a/Fe] ratios than the QSO-DLAs. The question that fol- 
lows is whether these higher values indicate enhanced a 
abundances (nucleosynthesis) and/or a higher depletion 
level for Fe (dust). 

From the nucleosynthetic viewpoint, one expects en- 
hanced a/Fe in the gas near GRBs because (i) the pro- 
genitors are massive stars and (ii) the high specific star 
formation rates of their host galaxies imply ages that are 
young compared to the tim e-scales for Type la enrich- 
ment (Christcnscn. Hiorth fc Gorosabel 2004). We con- 
tend that [a/Fe]> -f 0.3 dex is strongly expected for GRB- 
DLA gas from nucleosynthesis enrichment alon e (see also 
ICalura. Dessauges-Zavadskv fc Prochaskall2007f) . Regard- 
ing dust, one may expect high depletion levels for gas in 
or near star forming regions. This could also explain the 
offset in a/Fe between the GRB-DLAs and QSO-DLAs, 
especially if GRB-DLAs have systematically higher metal- 
licity (§[33). 

In the next subsection, we will consider dust in greater 
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Fig. 5. — Observed gas-phase relative abundances as a function of alpha abundance [a/H], zinc abundance [Zn/H], or metallicity [M/H]. 
The four panels show the values for GRB-DLAs and QSO-DLAs (sample HR-E) for (a) [a/Fe] ratios, (b) [Zn/Fe] ratios, (c) [Ti/Fe] ratios, 
and (d) [N/a] ratios. The GRB-DLAs are characterized by large [«/Fe] and [Zn/Fe] ratios indicating significant nucleosynthetic enhancement 
by massive stars and/or differential depletion. One also notes several solar and sub-solar upper limits on [Ti/Fe] which are strong evidence 
that refractory metals are depleted onto dust grains in GRB-DLAs. 



deptli. Before proceeding, we wish to emphasize the rela- 
tively low a/Fe value for GRB 050730: [a/Fc]=0.25 ± 0.15 
at [a/H] = —2.25. If we were to interpret this a/Fe ratio 
in terms of depletion, the intrinsic (nucleosynthetic) ratio 
would be approximately solar or even sub-solar; this would 
require star-form ation with low efficiency o ver time-sca les 



the level of a few tenths dex fe.g. iProchaska et al. ll20do( ). 
large Zn/Fe enhancements - especially at high metallicity 
- are unlikely to be explained by nucleosynthesis alone. 

We present the Zn/Fe ratios of GRB-DLAs and QSO- 
DLAs in FigurejSJa. Similar to the a/Fe ratios, the Zn/Fe 
values for the QSO-DLAs increase with metallicity. Again, 



of a f ew lOOlVIyr (jCalura. IVIatteucci fc Vladilo 2003 ; Dessauges-Zlai gaitebvisitfcfiA) -ained by higher depletion levels at higher 

|2007( ). This mode of star formati on does not appear com- metallicity. The GRB-DLAs do not exhibit this trend; the 

mon for high z GRB host galaxies (jChristensen. Hiorth fc Gorosabki es are uniformly large. This suggests that the majority 



[200I . In this one case, we argue that the gas is essentially 
undepleted. One may speculate further that the majority 
of QSO-DLAs at low metallicity are also nearly undepleted 
and that their observed a/Fe ratios simply imply signifi- 
cant a-enrichment. 

3.3.2. Differential Depletion 

It is routine in QSO-DLA studies to examine the Zn/Fe 
ratio to gauge the level of differential depl etion (or Zn/Cr; 
llVlever fc RothI Il990t [Pettini et al.\\l994 . As discussed 
above, this assumes that the nucleosynthetic production 
of Zn tracks the Fe peak closely such that enhancements 
in the gas- phase abundance of Zn/Fe is primarily due to Fe 
depletion. While this assumption should be questioned at 



of GRB-DLAs have large dust-to-metal ratios. Given the 
uncertainty on the nucleosynthesis of Zn, however, it is 
best to address this issue from yet another angle. 

In Figure [5j; we present upper limits to the Ti/Fe ratios 
of those GRB-DLAs with observations of the Till A1910 
transitions. In stellar atmospheres, Ti behaves like an 
a-element (i.e., Ti tracks Si,0,]Vtg in metal-poor stars; 
|]V[cWilliamlll997l : iProchaska et a/. II2OOOD . In the Galac- 
tic ISIVI, Ti is hi ghly refractory and one ge nerally observes 

[Ti/Fe]/sA/ < (|Jenkinslll987f ). Therefore.'Dessauges-Zavadsk v. Pr 



emphasized that observations of Ti/Fc lend to degenerate- 
free interpretation: super-solar Ti/Fe ratios indicate a- 
enhancement whereas sub-solar Ti/Fe ratios require differ- 
ential depletion. Therefore, the GRB-DLAs with [Ti/Fe] < 
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dex imply substantial depletion, especially in light of the 
large a/Fe ratios. We consider this definitive evidence that 
at least some GRB-DLAs are highly depleted. 

The results in Figures [5Jd,c argue that the GRB-DLAs 
have at least modest depletion levels. We contend that dif- 
ferential depletion contributes at least -1-0.3 dex to [a/Fe]g 
for most GRB-DLAs, i.e. at least 50% or more of the Fe 
is locked into dust grains. It remains an open question, 
therefore, whether the intrinsic a/Fe ratios of GRB-DLAs 
are enhanced relativ e to solar as expected for young, star- 
formi ng regions fe.g. lCalura. Dessauges-Zavadskv fc Prochaskal 
I2007D . While common practice is to examine the a/Zn ra- 
tio (again, using Zn as a proxy for Fe) , these results are es- 
pecially subject to the nucleosynthetic history of Zn. Nev- 
ertheless, examining a/Zn we find that the GRB-DLAs 
primarily show lower limits that are consistent with the so- 
lar abundance but allow for large enhancements. In pass- 
ing, we note that a large dataset of a/Zn measurements 
for the GRB-DLAs could inform the processes of Zn pro- 
duction and also interpretations of these measurements in 
QSO-DLAs (|Nissen et al. l[200l . 

3.3.3. N/a 

Another abundance ratio of particular interest in terms 
of nucleosynthesis is N/a. Because N is believed to be pro- 
duced primarily in the AGB phenomenon of intermediate 
mass stars (e.g. iMevnet fc Maedeil I2002D . the N/a ratio 
provides a diagnostic of the sta r formation history of the 
galax y, especially at early times (|Henrv. Edmunds fc KoppenI 
|2000( ). Atomic nitrogen is the dominant ionization state 
in neutral regions and, unfortunately, this low-ion's transi- 
tions all lie within the Lya forest. Therefore, the N abun- 
dance is difficult to measure in absorption-line systems: 
high resolution data is required and line-profile analysis is 
often required to recover N° column densities. Because N 
is non-refractory, the observed ratios should nearly corre- 
spond to the intrinsic values. The only s erious systematic 
uncer tainties are ionization corrections (|Prochaska et al.\ 
I2002D which should be small for GRB-DLAs given the very 
large H I column densities. 

The current GRB-DLA sample only allows one to con- 
strain the NO column density for three GRB-DLAs: GRB050730 
GRB050820, GRB050922C. The N/a ratios are presented 
in Figure [SJi in comparison with measurements made for 
the QSO-D LAs (Prochaska et al. "200j; ICenturion et al. I 
120031 : iDessa ugcs-Zavadskv et al. 200^. Although this is 
a small sample, it is evident that there is a large dispersion 
in [N/a] for the GRB-DLAs. This indicates that their host 
galaxies have experienced a diverse range of star formation 
histories. 

Two of the N/a values (GRB 050820, GRB 050730) bracket 
the locus of observations for the QSO-DLAs. These values 
are relatively high ([N/a] > —1) indicating a significant 
enhancement by intermediate mass stars. This is not sur- 
prising for GRB 050820 where the gas metallicity is large, 
implying several generations of star formation. The N/a 
value for GRB 050730 (at met allicity [M/H] » -2) sug- 
gests an age older than 200 Myr (jHenrv. Edmunds fc KoppenI 
l2000t ). which at z = 4 implies a formation redshift ap- 
proaching the epoch of rcionization. 

The most intriguing measurement, however, is the up- 
per hmit to N/a in GRB 050922C. This limit lies well 



belo w even the so-called lower plateau of QSO-DLA val- 
ues (jProchaska et al. 1 120021 : ICenturion et al. I l2003f) . If 
this upper limit has a nucleosynthetic origin, it implies 
gas that has not been polluted by any intermediate mass 
stars. Given this rather extreme value, it is important 
to further address ionization corrections; perhaps the low 
N/a value indicates an unusual ionization state for the gas 
(i.e. a high N^/N" ratio). Because ato mic nitrogen has 
a rela tively large cross-section to X-rays (|Sofia fc JenkinsI 
IT998I) . it may be under-abundant even in neutral regions. 



Both the GRB afterglow and the progenitor imply an en- 
hanced radiation field which could lead to high N^/N^ 
ratios. We are especially concerned about photoioniza- 
tion for GRB 050922C given that it has a low MgO col- 
umn density and likely a, bsorption by the N V doublet 
(jPiranomonte et al. I l2007l ) . Equilibrium photoionization 
calculations, however, indicate ionization corrections are 
generally less tha n 0.5 dex, even for an in put spectrum as 
hard as a quasar (|Prochaska et al. II2002D . In a future pa- 
per, we will perform a time-dependent calculation to track 
the effects of the afterglow on this and other ionization 
states relevant to our observations. 

Unless the ionization corrections are larger than 0.5 dex, 
the limit on [N/a] lies below that of any other astrophys- 
ical environment. This is not especially surprising for a 
GRB event, however, especially one with very low metal- 
licity. In this case, the galaxy may have initiated star 
formation only very recently and the observed metals may 
have no contribution from intermediate mass stars. For a 
standard initial mass funct ion, this requires an enrich ment 
age of less than - 30 Myr (jHenrv fc Prochaskall2007[ ) . Fu- 
ture observations will reveal whether low N/a values are 
characteristic of a larger sample of GRB-DLAs. 
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Fig. 6.— C°/C+ ratios of the GRB-DLAs (dark) and QSO-DLAs 
(light) inferred from column density measurements of C I transitions 
and by assuming logA''(C+) = [a/H] + log A^jjl — 4. Typical values 
of the cold n eutral medium of the Galactic ISM are 0^/0+ = 10"-^ 
llJenkins fc S hava 1979) while significantly lower values (10"'*) cor- 
respond to the warm neutral medium (Liszt 200^). We interpret 
the observations as more reflective of the latter phase of the ISM, 
but caution that the results are sensitive to the far-UV intensity, 
metallicity, and cosmic-ray flux through the gas. 
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3.4. Atomic Carbon 

In the Galaxy, an excellent tracer of cold, dense gas is 
atomic carbon. Because C° has an ionization potential 
below 1 Ryd, it is generally photoionized by the ambient 
far-UV radiation field of the ISM. In cold and dense re- 
gions, however, the recombination rate (proportional to 
the electron and carbon densities) is sufficiently large that 
cold clouds frequentl y show detectable abs orption from 
C I transitions fe.g. iJenkins fc Tripd [200ll ). For stan- 
dard Galactic ISM conditions, the C^/C"*" ratio is w 10"^ 
in the cold neutral mediu m and decre ases to « 10^'' in 
the warm neutral medium ([Lisztj [2003 ) . These ratios, of 
course, are sensitive to the gas metallicity, f ar-UV radi- 
ation field, and cosmic r ay flux of the ISM (|Lis^ [200l 
IWolfe. Prochaska fc Gawiser 2003 ) . Nevertheless, very low 
C°/C"'" values generally imply gas with density less than 
10 particles per cm'^. 

Figure [6] presents the C"/C+ ratios estimated for the 
GRB-DLAs and QSO-DLAs (sample HR-A). While we can 
measure NlC'^) directly from C I transitions, the C II 
transitions are generally too saturated to provide even 
a valuable lower limit to iV(C+). Therefore, we infer 
A^(C+) from the [a/H] abundance and iVni v alue, assum- 
ing [q ;/C]=+0.3 (i.e. a modest a-enhancement: lAkerman et al. 
[200I : 

log N{G+) = [a/H] + log Nm - 12 + 8.3 . (2) 

Note that this calculation is actually independent of the 
A^Hi value. None of the GRB- DLAs in our sampl e has a 
positive detection, although Vreeswiik et al. I ()2006l ) report 
the detection of C I in GRB-DLA 060418 with a value con- 
sistent with our upper limit. Because the metal-line col- 
umn densities of the GRB-DLAs are large, the upper lim- 
its to C°/C+ generally lie below those of the QSO-DLAs 
even though the latter are derived from higher S/N data. 
The observations suggest the gas in GRB-DLAs (and many 
QSO-DLAs) is dominated by a warm, less dense phase. We 
caution again, however, that the predicted G^/C~^ values 
are sensitive to the far-UV flux, which is likely to be en- 
hanced near GRB. A more cautious interpretation of the 
data is that it offers additional, yet non-conclusive support 
that GRB-DLA gas is not characteristic of the Galactic 
cold neutral medium. 

4. DISCUSSION 

4.1. General Characteristics of the ISM Near Star 
Forming Regions 

Reviewing the results presented in the previous section, 
wc summarize the properties of gas observed surrounding 
GRB afterglows as having: 

• Very large H I column densities with median 
Nm ~ 1021-6 cm-2. 

• A wide range of metallicities (1/100 to nearly 
solar) with median larger than 1/10 solar. 

• Observed a/Fe ratios in excess of 3x the solar 
abundance reflecting enrichment by massive stars 
(Type II SN) and/or differential depletion. 



• Large [Zn/Fe] values and several solar or sub-solar 
[Ti/Fe] upper limits which indicate substantial 
differential depletion. 

• A dispersion of [N/q] indicating a diverse set of 
star formation histories. We identify one case 
(GRB 050922C) with an extremely low value 
suggesting gas enriched solely by massive stars. 

At first glance, the observation of large H I surface 
densities and significantly depleted, chemically enriched 
gas appear to reflect the conditions expected for actively 
star- forming regions (i.e. molecular clouds). Quantita- 
tively, however, we find that the GRB-DLA properties are 
more characteristic of the ambient ISM of modern galaxies. 
As several authors have noted for individua l GRB-DLA 
(ISavagho fc Falll l2004l : iPenprase et~ar\ 120061: IShin et al. I 
I2006D . the observed depletion levels are more representa- 
tive of w arm (i.e. less d ense) clouds in the Galactic ISM 
(see also ISavagliol l2006l) . On the other hand, the major- 
ity of GRB-DLAs have sub-solar metallicity and a direct 
comparison with the Galactic ISM may not be appropri- 
ate. Unfortunately, the differential depletion characteris- 
tics of cold, dense gas in other galaxies (e.g. the Magel- 
lanic clouds) has not been extensively measured, in part 
because dust obscuration challenges the analysis. In the 
Small Magellanic Cloud (SMC), one does observe highly 
depleted gas along the sightline to Sk 155 but one also finds 
modest depletions in clou ds where C I analysis indicates 
large density (> 100 cm '^ IWeltv et aLll200ll : ISofia et al. I 
l2006tlWeltv al. Il2006t ). We cautiously conclude that the 
depletion levels observed for the GRB-DLAs are generally 
lower than that observed for cold, dense clouds in the lo- 
cal universe. It would be very valuable to independently 
estimate the density of gas responsible for GRB-DLAs. 

Throughout this paper, we have presented additional 
observations which argue the GRB-DLAs do not arise in 
a cold, dens e phase. First, the gas has very low molecu- 
lar fraction (jVreeswiik et al. |[200l. Paper III). Although 
a low molecular fraction could result from photoionization 
by the progenitor, nearby O and B stars, and the GRB 
afterglow, we have argued the gas lies at > 100 pc from 
the GRB where these effects may be lessened. If the gas 
has large particle density {uh > lO^cm"'^), then it would 
be difhcult to maintain such a low H2 fraction (Paper III). 

Second, the only atomic line routinely measured in the 
GRB-DLAs is Mg I and its large dielectronic recombina- 
tion coefficient does not require cold, dense gas. In ji l3.4l we 
presented upper limits on C^/C"*" which argue against gas 
with characteristics like the dense CNM of the Milky Way. 
Third, there are no especially anomalous abundance pat- 
terns observed for the gas (e.g. extreme metallicity, very 
high a/Fe ratios) which could occur if the gas represented 
partially mixed SN ejecta or circumstellar material. Al- 
together, we conclude that the gas revealed by GRB af- 
terglow spectra represents the ISM near the GRB but not 
gas directly associated with its own molecular cloud. 

4.2. Comparisons with QSO-DLAs 

In the previous subsection, we argued that the GRB- 
DLAs represent gas from the ambient ISM of their host 
galaxies. The large A^hi values of QSO-DLAs argue these 
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sightlines also penetrate the ISM of high z galaxies. Fur- 
thermore, the majority of QSO-DLAs have depletion pat- 
terns, low C+O/C"'" ratios, and low molecular fractions 
which are not characteristic of cold, dense gas in the Galaxy. 
In this respect, it is fair to compare the observations of 
the two populations to address whether they are drawn 
from the same parent population of galaxies. In § [3] wc 
compared the A^hi values, metallicity, and relative abun- 
dance ratios of the QSO-DLAs and GRB-DLAs. The lat- 
ter are characterized by higher A^hi values, higher a/Fe 
and Zn/Fe ratios, and higher metallicities. These differ- 
ences beg the question: Are GRB-DLAs and QSO-DLAs 
drawn from distinct distributions of host galaxies ? 

We contend that the observations presented in this pa- 
per do not contradict the null hypothesis that GRB-DLAs 
and QSO-DLAs have identical parent populations of galax- 
ies but their differences can be explained by sightlines with 
distinct impact parameter distributions. Although offsets 
exist in [M/H] and the X/Fe ratios, the differences are 
relatively small (< 0.5 dex). Typical observed gradients 
in the metallicity of stars and H II regions of local galax- 
ies imply a variation of at least 0.3 dex from the center of 
the galaxy to the r adii roughly corresponding to Ahi = 
2 X 10^° cm"^ (e.g. iKennic utt. Bresolin fc GarnettI 120031: 



l2005f ). This is most easily seen by extrapolating the /(A^hi) cx 
N^^ power-law observed at A^hi < 10^^ cm~^ to larger 
column densities. One predicts that there should be « 10 
damped Lya systems with A'hi > 10^^ cm~^ per 500 ob- 
served. The absence of a single QSO-DLAs with this 
column density indicates that /(A^hi) has a steeper than 
A'jjj^ dependence at these column densities. The observa- 
tion of a break in /(A^hi) and the routine observation of 
A^Hi > 10^^ cm~^ in GRB-DLAs raises concerns that dust 
obscuration is biasing the QSO-DLA sample against large 
A^Hi values. 

There is additional (circumstantial) evidence for a dust 
obscuration bias in QSO-DLAs. First, QSO-DLAs with 
a combination of large A^hi and high metallicity (metal- 
strong DLA), i.e. potentially large dust column s, are rare 
(jBoisse et aL 1 119981 : [lerbert-Fort et a/. Il2006t) . Second, 
every DLA i s enriched by heavy m etals with metallicity 
[M/H] > -3 (|Prochaska et al. II2003D and many have rela- 
tive abundance rat ios (e.g. Z n/Fe) indicative of differential 
dust depletion (Pet tini et al. 1994; Prochaska & Wolfe 200|). 
Third, several quasars with intervening DLA show signif- 
icant reddening that can be attributed to the absorber 
IVladilo et al. I ()2006D . Motivated by these observations, 
we will test whether QSO-DLA samples would be likely to 



Chen. Kennicutt fc Rauch 200i|Luck, Kovtyukh fc Andrievskytndude GRB-DLA sightlines if background quasars were 
2006D. Furthermore, the H I surface density and volume placed directly on the sightline 



density are undoubtedly larger toward the center of the 
galaxy (at least until the gas becomes molecular). And it 
is reasonable to expect that higher metallicity, denser gas 
will exhibit larger depletion levels than more tenuous gas 
at the outer edges of H I disks. 

At present, the only distribution with an especially large 
offset (> Idex difference in median value) is the A'hi val- 
ues. While differences in the impact parameters of GRB- 
DLAs and QSO-DLAs may explain the offsets in the me- 
dian of the distributions for A'hi, [M/H], and X/Fe ra- 
tios, the identification of many GRB-DLAs with A^hi > 
10^^ cm~^ and none in the QSO-DLA sample is striking. 
In § 13. 1[ we demonstrated that the GRB-DLA A^hi dis- 
tribution is reasonably well explained by an H II region 
embedded with an exponential H I disk with central col- 
umn density logA^o = lO^^cm"^. While this simplistic 
model is a good description of the observations, it does lead 
to a conflict with the H I frequency distribution /(A'hi) 
observed for the QSO-DLAs. If the majority of high z 
galaxies have logA^o = lO^'^cm^^, then one would not 
observe a break in /(A^hi) at A^hi ^ lO^^'^cm "^ as ob- 
served (iProchaska. Herbert-Fort fc Wolfdl2005l ) . One pos- 
sible resolution for this conflict is that GRB-DLA sight- 
lines are obscured by dust from magnitude-limited QSO 
surveys. Let us now consider this hypothesis in greater 
detail. 

4.3. Can QSO-DLA Samples Probe GRB-DLA 
Sightlines ? 

As Figure[2]reveals, a significant fraction of GRB-DLAs 
have A^Hi > 10^^ cm'^ yet not one of the ^ 500 QSO- 
DLAs discovered thus far have A^hi > lO^^'^cm"^. An 
analysis of the SDSS dataset actually requires that the 
QSO-DLAs exhibit a break in their H I frequency distribu- 
tion /(A^hi) at A'hi ~ 1 021-5 Qi^-'^ explain the absence 
of large A^hi absorbers (jProchaska. Herbert-Fort fc Wolfd 



We begin by deriving the dust-to-gas ratios of the GRB- 
DLAs. In the following, we will assume the GRB-DLA 
dust properties are best matched by those observed for 
the SMC. We make this choice because (i) no GRB-DLA 
to date has exhibited the 2 175 A dust feature character- 
istic of the M ilky Way a nd LMC extinction la ws (e.g. 
| Mira.bal et 011,2002: Sa va gho fc Fal l 2004; EUiso n et aU 
' 2006 1: iButler et al. Il2007t ): and (ii) the metaUicities of GRB- 
DLAs are representative of the SMC. We derive the dust- 
to-gas ratio of the DLAs relative to the SMC: k/ksmc 
by assuming the SMC has metaUicity [M/H]5a/c = —0.7 
and that 90% of its refractory element s are depleted from 
the g as-phase, i.e. [M/Fe]sA/c = +1 (|Weltv et al. Ill997l 
l2001f ). In this case, the relative dust-to-gas ratio is 

- 1 _ 10(^'-[WFe]) - 



KSMC 



1 - 10- 



(3) 



where allows for a nucleosynthetic contribution to the 
observed [M/Fe] value and we assume = for the SMC. 
For the GRB-DLAs, we adopt the [M/H] and [M/Fe] val- 
ues listed in Tabled In the majority of cases, M is sulfur 
and the remainder assume Zn or Si. In the following, we 
will assume A^ = 0.3 dex for the GRB-DLAs and 0.2 dex 
for the QSO-DLAs. The high specific star formation rates 
observed for GRB host galaxies indicate young ages and, 
therefore, gas enriched by massive stars giving a-enriched 
abundances ([a/Fe] > -F0.2). Similarly, the QSO-DLAs 
show a plateau of [Si/Fe] =-1-0.25 dex values at low metal- 
licity indicative of Type II enrichment (jProchaska fc Wolfd 
2002 ), and enhanced a/Fe ratios in "dust free" systems 
(jDessauges-Zavadskv et aL"l 120061 ). In any case, the re- 
sults are not sensitive to our choices for A^ unless we allow 
A, > 0.5 dex for the GRB-DLAs or A^ > 0.3 dex for the 
QSO-DLAs. 

Figure [7] shows the k/ksmc values for the GRB-DLAs 
and also the echelle sample (HR-E) of QSO-DLAs. Both 
samples exhibit a wide range of dust-to-gas ratios relative 



14 



10.000 



1.000 



J 0.100 
■2 



0.010 



0.001 





. GRB-DLA ; 




, QSO-DLA ■ 


- + + 








^ t : , , ' * 








; + + * 









20.5 21.0 21.5 22.0 22.5 23.0 
logNm 



0.20 



0.15 



u 

S 0.10 



0.05 



0.00 



-4 ,tt't 



Hf4++ +f + M 



20.5 21.0 



21.5 



22.0 22.5 



Fig. 7. — Dust-to-gas ratio k of the GRB-DLAs (dark blue) and 
QSO-DLAs (light green) relative to the dust-to-gas ratio of the SMC 
i^SMc) a function of H I column density. 



to the SMC. We also find that the GRB-DLAs show a 
higher fraction of large values {k/ksmc ^ !)■ Using the 
observed TVhi values and assuming SMC dust properties, 
we can convert these dust-to-gas ratios into visual extinc- 
tion, Av- Specifically, we adopt two e mpirical relations: 
(i) th e observed reddening in the SMC (|Tumlinson et a/.~l 
120021 ) 

log E{B -V)^ log Nm ~ 22.95 (4) 



and ( ii 
I200I ] 



) a linear relation between Ay and E{B—V) ([Gordon et 
modified by the dust-to-gas ratio 



Fig. 8. — Predicted visual extinctions Ay for the GRB- 
DLA and QSO-DLA sightlines assuming the dust-to-gas ratios pre- 
sented in Fi gure [7] and assumi ng the dust pr operties of the SMC 
IITumlinson et al. 2002; Gord on et al. I [200^ . All of the QSO- 
DLA sightlines show very small values {Ay < 0.05 mag). In con- 
trast, the GRB-DLA show a nearly bimodal population with several 
small values and several sightlines approaching 0.2 mag. Note that 
Ay = 0.1 mag corresponds to an extinction of 0.6 mag assuming the 
SMC extinction law. 



against their detection in QSO samples. To estimate the 
effects, we assessed the signal-to- noise ratio (S/N) at the 
center of the Lya profile of each QSO-DLAs in the SDSS- 
__DlR3 sample. Specifically, we measured the median S/N 
-Selli ng the continuum adopted in the fi t s to th e Lya profiles 



Ay = 2.74 E{B - V) 



(5) 



KSMC 

FigurelUpresents the Ay values against H I column den- 
sity for the GRB-DLAs compared against QSO -DLAs. All 

of the QSO-DLAs have Ay < 0.05 mag (see also lProchaska fc Wolfd 
I2OO2I) . It is not surprising, therefore, that one obs erves 
very little reddening of quasar spectra by QSO-DLAs (jMurphv 
l2004f ). In contrast, the estimated Ay values for approxi- 



bv lProchaska. Herbert-Fort fc Wolfd (|2005f ). By the defi- 
nition of their statistical sample, the S/N exceeds 4pix~-'^ 
in each case. We then lowered the S/N by adopting an ex- 
tinction Ay assuming that the noise is quasar dominated, 



S/N= (S/No)/VlO^^«°/2-5 (6) 
The fra ction of QSO-DLAs that satisfy the S/N> 4pix~i 



mately half the GRB-DLAs are substantially larger than 
the QSO-DLA distribution {Ay > 0.1 mag). At Ay 
0.1 mag, the extinction at the Lya transiti on is A^yg — 
0.6 m ag assuming the SMC extinction law (jPrevot et al. I 
Il984[ ). This is a considerable but not overwhelming level 
of extinction. Interestingly, there is no obvious correla- 
tion between Ay and iVni, at least for Ay values greater 
than 0.05 mag. We also note that the Ay values presented 
in Figure [5] are roughly consistent with the Ay values in- 
ferred from broad-band photometry or spec trophotometry 



of GR B afterglows (e.g.lM irabal et al. 2002| ; | Vreeswiik et al. I 
[2OOI ISavaglio fc Fai]n20b4l: iButler et al. I I2006D . On "the 
other hand, our Ay values appear to be at odds with 
those inferred from X-ray absorption m easurements (e.g. 
IWatson et al. Il2006l : l&ltler aZ. II2006D which has led some 
authors to invok e "grey" extinction laws for the dust en- 
veloping GRBs (jCalama fc Wiiersll200lD . This issue will 
be addressed in greater detail in a future paper (Butler et 
al. 2007, in prep). 

Finally, we examine whether the inferred Ay values 
for the GRB-DLA sightlines are sufficiently large to bias 



fc^jjtskiBjn is shown in Figure [9] as a function of Ay. For 
Ay < 0.2 mag, dust obscuration has a relatively small im- 
pact on the QSO-DLA sample. We cautiously conclude, 
therefore, that the majority of observed GRB-DLAs would 
also be observed along quasar sightlines. 

This conclusion, however, is subject to a few caveats. 
First, the exercise described by Figure [9] does not account 
for quasars that would become too red or too faint to sat- 
isfy t he QSO target criteria for the SDSS (jRichards et al. I 
' 2004t ). We suspect, however, that this effect is much smaller 
than the S/N cut we have applied (Murphy et al. 2007, in 
prep). Second, several of the k/ksmc values for the GRB- 
DLAs are lower limits because of line saturation. It is pos- 
sible that these have Ay > 0.4 mag and would have a small 
probability of entering the QSO-DLA sample. Third, all 
of the results assume SMC extinction. This assumption 
has two implications: (i) it establishes the reddening per 
H I atom (scaled by metallicity; Equation [3]) and (ii) the 
shape of the extinction law which determines A^ya/Ay. 
Regarding the latter aspect, the SMC has the steepest ex- 
tinction law of local galaxies. Other ("grayer") dust laws 
would generally lead to less extinction at Lya. The red- 
dening per H I atom in the SMC, however, is smaller than 
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Av [SMC] 

Fig. 9. — Fraction of DLAs recovered from the SDSS QSO-DLA 
survey of IProchaska. He rbert-Fort & Wolfe (2005) as a function of 
additional visual extinction (Ay) for the DLAs. For Ay < 0.2 mag, 
the effect of obscuration is mild while only a small fraction of QSO- 
DLAs with Av > 0-4 mag would be recovered in the SDSS survey. 



the Milky Way and not only b ecause of the SMC's lower 
metallicity (|Weltv et al. Il2006f ). Fourth, we have ignored 
the possibility that a significant column of dust exists in 
the H II region surrounding the GRB. This could be re- 
vealed by broad-band photometry of the afterglow com- 
pared against Ay estimations from the observed metal 
column densities. Finally, we note that a non- negligible 
fraction of GRB afterglows are sever ely reddened, most 
likely by dust surro unding the event (|Levan et al. Il2006t 
iPellizza et al. I [200^ . These "dark bursts" do not enter 
our GRB-DLA sample and likely have Ay values much 
larger than those presented in Figure [8] At present, this 
extremely dusty gas is not well probed by GRB or QSO 
observations. 

These concerns aside, our analysis suggests that many 
(if not all) of the GRB-DLAs in our sample could be 
detected in modern QSO-DLA surveys. If confirmed by 
larger samples, then the very low incidence of QSO-DLAs 
with A^Hi > 10^^ cm~^ must be explained by a substan- 
tial decrease in the cross-section of gas with these column 
densities. To avoid conflict with the observed /(iVni) dis- 
tribution, one requires that the cross-section of the re- 
gions probed by GRB-DLAs must be less than 1% of the 
cross-section corresponding to the damped Lya criterion, 
log A^Hi = 20.3. Characterizing the latter by a radius r2o.3, 
we require the vcrb < '^20.3/10- This relation is rather 
easily achieved if rcRB is lOOpc or less, but cou ld pose a 
problem if rcuB ^ Ikpc (jVreeswiik et al. |[2006D . 

We wish to thank Arthur M. Wolfe for helpful discus- 
sions and for allowing us to analyze results on QSO-DLAs 
prior to the public release of these data. We thank E. 
Ramirez-Ruiz for constructive comments on an early ver- 
sion of this manuscript. J.X.P. is partially supported by 
NASA/Swift grant NNG05GF55G and an NSF CAREER 
grant (AST-0548180). 
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